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ABSTRACT 


Scaling  of  noise  generation  due  to  a  single  bubble  in  a  cavitating  trailing 
vortex  is  studied.  Vortex  core  dynamics  between  full-scale  and  model  is 
mathematically  formulated  based  on  a  similarity  flow  approach  and 
McCormick’s  hypothesis.  Acoustic  frequency  and  amplitude  scaling  formulas 
to  obtain  full-scale  acoustic  spectrum  from  model  data  are  derived.  The  acoustic 
scaling  formulas  are  also  compared  with  bubble  dynamic  calculations  with 
encouraging  results.  The  numerical  calculations  are  obtained  from  a  Dynaflow 
Inc.  computer  code.  The  acoustic  scaling  formulas  are  compared  with 
Strasberg’s  classic  X  scaling  and  noise  level  scaling. 


ADMINISTRATIVE  INFORMATION 

The  work  was  partially  supported  by  NAVSEA  93R  under  the  direction  of  Doug  Dahmer 
and  Jude  Brown  and  partially  supported  by  NAVSEA  PMS  450  Virginia  Class  Propulsor 
Program  under  the  direction  of  Larry  Becker  and  Jack  Lee.  Work  unit  numbers  are  02-1-5080- 
218-52  and  02-1-5080-228-69. 


INTRODUCTION 

A  three-dimensional  lifting  surface  has  a  vortex  trailing  from  its  tip  caused  by  the 
crossover  of  fluid  from  the  high-pressure  side  to  the  low  pressure  side.  For  propellers,  the 
vortices  from  the  blades  form  a  helical  path  downstream  of  the  blades.  The  vortices  spin  faster 
and  cause  regions  of  low  pressure  in  their  core  at  conditions  of  higher  speeds  or  higher  propeller 
thrust  (loading).  When  the  pressure  in  the  core  drops  below  the  vapor  pressure  of  the  flowing 
liquid,  cavitation  can  occur.  On  a  marine  propeller,  cavitation  typically  occurs  in  the  trailing 
vortex  before  it  occurs  on  the  blade  surface.  The  prediction  of  the  onset  of  tip  vortex  cavitation 
is  of  major  interest  to  the  Navy  because  of  the  noise  and  hull  vibration  associated  with  the  vortex 
cavitation.  The  flow  field  in  a  trailing  vortex  is  sufficient  complex  that  prediction  of  cavitation 
inception  and  noise  generation  still  relies  heavily  on  model  tests  extrapolated  to  full-scale  using 
cavitation  scaling  laws. 

When  a  model  propeller  is  tested,  cavitation  inception  numbers  Oj  with  respect  to 
different  ship  speeds  and  depths  are  determined.  Acoustic  signatures  associated  with  various 
stages  of  cavitation  are  measured.  Three  scaling  formulas  are  needed  to  translate  model  data  to 
full-scale  performance;  (1)  cavitation  inception  scale  to  predict  cavitation  inception  speed  of  the 
full-scale  propeller  at  a  given  depth  or  cavitation  inception  depth  at  a  given  speed,  (2)  frequency 
scaling  of  acoustic  spectra,  and  (3)  amplitude  scaling  of  acoustic  spectra.  In  this  report, 
frequency  scaling  and  amplitude  scaling  of  noise  generated  by  bubble  cavitation  are  studied. 
The  cavitation  inception  scaling  has  been  formulated  and  presented  by  Shen  et  al  [2001],  Shen, 
Jessup  and  Gowing,  [2001],  Hsiao,  Chahine  and  Liu  [2000]  and  Hsiao  and  Chahine  [2002]. 


'  References  are  listed  in  alphabetical  order  on  page  25. 
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Cavitation  noise  and  associated  noise  scaling  has  been  studied  by  many  researchers  (see 
review  papers  by  Strasberg  [1977],  Chahine  [1979],  Blake  and  Sevik  [1982],  Blake  [1986]  and 
Baiter  [1982,  1989]).  Due  to  complexity  of  the  flow  field  in  a  tip  vortex,  the  effect  of  pressure 
and  vortex  core  dynamics  on  bubble  cavitation  and  noise  generation  has  not  been  considered. 
Experimental  data  indicate  that  tip  vortex  cavitation  inception  is  sensitive  to  the  Reynolds  scale 
(McCormick  [1954])  and  free  stream  nuclei  (Acosta  and  Parkin  [1975]).  Scaling  of  noise 
generation  due  to  bubble  cavitation  in  trailing  vortices  is  studied  in  this  report. 

In  the  first  part  of  this  report,  vortex  core  dynamics  between  full-scale  and  model  is 
mathematically  formulated.  The  scaling  formulas  to  obtain  full-scale  acoustic  spectrum  form 
model  data  are  derived.  The  present  acoustic-scaling  formulas  are  compared  with  Strasberg’s 
classic  X  scaling  and  noise  level  scaling.  In  the  second  part  of  this  report,  the  present  acoustic¬ 
scaling  formulas  are  evaluated  with  bubble  dynamic  calculations.  The  numerical  calculations 
are  based  on  the  modified  Surface  Averaged  Pressure  (SAP)  Rayleigh-Plesset  bubble  dynamic 
equation  (Hsiao  et  al  [2002]).  Bubble  growth,  collapse,  and  radiated  noise  due  to  a  cavitation 
bubble  are  computed  for  a  trailing  vortex  with  prescribed  circulation  and  pressure  distributions 
along  the  vortex  core. 


PROBLEM  FORMULATION  AND  SOLUTION  APPROACH 

To  model  the  vortex  and  roll-up  process,  consider  a  propeller  blade  with  a  system  of 
trailing  vortices  as  shown  in  Figure  1.  Let  V®,  and  P®  denote  the  free  stream  velocity  and 
reference  pressure,  y  the  local  bound  vorticity  on  the  foil.  Let  Rp  and  c  denote  the  propeller 
radius  and  local  chord  length.  Let  T  represent  the  vorticity  in  the  trailing  vortex,  Cp  is  the 
pressure  coefficient  in  the  vortex  core,  and  ac  is  the  radius  of  the  vortex  core. 

Let  X  denote  the  linear-scale  geometric  ratio  of  full-scale  to  model-scale,  namely  Rpf 
=?JRpm  and  Cf  =  X  cm,  where  subscripts  ‘m’  and  ‘f  denote  the  model  and  full-scale,  respectively. 
Let  J  denote  the  propeller  advance  coefficient.  In  the  scaling  study,  it  is  assumed  that  the 
advance  coefficients  are  the  same  for  full-scale  propeller  and  model  propeller.  The 
hydrodynamic  angle  into  the  blade  section  will  be  the  same  for  model  and  full-scale.  Assuming 
both  model  and  full-scale  have  the  same  planform,  camber  and  section  profiles,  geometric 
similarity  between  the  model  and  full-scale  is  preserved.  For  uniform  inflow  and  same  angles  of 
attack  on  both  propeller  blade  section,  namely  <Xf  =  am,  kinematic  similarity  is  also  satisfied.  For 
dynamic  similarity,  it  is  assumed  that  the  trailing  edge  vortex  roll-up  process  between  model  and 
full-scale  is  similar.  For  simplicity  a  Rankine  vortex  model  will  be  used  to  relate  circulation, 
core  dimension,  tangential  velocity  and  pressure  distributions. 

The  process  of  bubbles’  capture  into  the  tip  vortex  is  complex.  In  the  present  study  it  is 
assumed  that  a  bubble  enters  the  vortex  core  and  convects  downstream.  The  vortex  core  is  a 
region  of  low  pressure.  For  sufficiently  low  pressures,  a  bubble  entering  a  vortex  core  will 
cavitate,  experience  growth  and  collapse,  and  radiate  acoustically.  It  is  presently  assumed  that 
the  bubble  remains  spherical  throughout  the  cavitation  process.  The  cavitation  spectrum  covers 
a  wide  range  of  frequencies.  The  scaling  of  peak  acoustic  frequency  and  amplitude  due  to 
collapse  of  a  single  bubble  in  a  trailing  vortex  is  investigated  here. 
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FREQUENCY  AT  PEAK  AMPLITUDE  DUE  TO  COLLAPSE  OF  A  SPHERICAL  BUBBLE 

The  behavior  of  a  spherical  bubble  in  a  pressure  field  for  incompressible  flow  is  given  by 
the  Rayleigh  equation  (Plesset  [1948],  Fitzpatrick  and  Strasberg  [1956]) 

R  (d2R/dt2)  +  3/2  (dR/dt)2  =  1/p  (pv  +pg-Pe  -2S/R)  s  1/p  (p(R)- Pe ),  (1 ) 

where  p(R)  =  pv  +pg  -2S/R  is  the  pressure  in  the  liquid  just  outside  the  bubble  surface.  R  is  the 
time  dependent  bubble  radius,  p  is  the  liquid  density,  pv  is  the  vapor  pressure  of  the  fluid,  pg  is 
the  gas  pressure  inside  the  bubble,  Pe  is  the  ambient  pressure  inside  the  vortex  core  local  to  the 
bubble,  and  S  is  the  surface  tension  of  the  fluid.  Under  tension  the  bubble  will  grow  as  it  is 
convected  downstream  inside  the  vortex  core.  In  a  high-pressure  zone,  the  bubble  will  collapse. 
Before  the  collapse,  the  bubble  attains  a  maximum  bubble  size,  which  is  denoted  by  Rmax-  It  is 
noted  that  p(R)  -  Pe  >  0  when  the  bubble  is  in  growing  stage  and,  p(R)  -  Pe  <  0  when  the  bubble 
is  in  collapsing  stage.  This  fact  is  utilized  in  the  following  derivation. 

For  the  collapse  phase,  Rayleigh  found  the  total  time  of  collapse  Tc  to  be  (Lamb  [1945]) 


T,  =0.9.5^,^.^))  p) 

where  Pe  -  p(R)  are  taken  when  R  =  Rmax.  Equation  (2)  is  only  approximate,  since  it  does  not  take 
into  account  the  fact  that  Pe  is  not  always  constant.  Ross  [1976]  presented  an  idealized 
cavitation  spectrum  as  shown  in  Figure  2  (also  see  Fitzpatrick  and  Strasberg  [1956]).  Measured 
cavitation  spectra  from  analyzing  the  available  experimental  data  show  the  peak  amplitude 
occurring  at  a  frequency  related  to  the  collapse  time  of  the  largest  bubble; 


fp=  !/2(l/Rmax  )/P(R)'P^  ,  (3) 

aside  from  difference  of  constants,  the  frequency  at  peak  amplitude  in  equation  (3)  gives  the 
same  expression  as  the  Rayleigh  bubble  collapse  time  in  equation  (2)  since  a  frequency  is  the 
inverse  of  time. 

Let  fpm  denotes  the  frequency  at  peak  amplitude  of  the  noise  spectra  generated  by  the 
model  due  to  a  bubble  cavitation.  Let  the  subscript  m  denote  the  variables  associated  with  the 
model.  We  have: 


pm 


(P(R)-Pe)„ 


(4) 


The  Rayleigh  bubble  dynamic  equation  expressed  in  Equation  (1)  is  quite  general  and  is 
applicable  model  scale  as  well  as  full-scale.  Let  fPf  denote  the  peak  frequency  of  the  noise 


3 


spectra  generated  by  the  full-scale  due  to  a  bubble  cavitation.  Let  the  subscript  f  denote  the 
variables  associated  with  the  full-scale.  We  have: 


fpf=’/*(l/Rma,r 


}J(P(R)-pcy/ 


From  (4)  and  (5),  we  have: 


V 


pm 


( Rmaxm  Y  V(P(R)-P,)r  Y 

AV(p(R)-pc)J, 


-Jpf 


(5) 


(6) 


The  frequency  of  the  noise  spectrum  generated  in  full-scale  due  to  collapse  of  a  cavitation 
spherical  bubble  can  be  predicted  from  the  model  data  by  Equation  (6).  Note  that  the  right  hand 
side  of  Equation  (6)  is  expressed  in  dynamic  variables.  To  be  useful  for  scaling,  the  dynamic 
variables  must  be  expressed  in  terms  of  geometric  and  kinematic  variables. 

The  Ratio  of  Rmax  m  /  Rmax /  in  a  Trailing  Vortex :  Experimental  data  indicate  that  tip 
vortex  cavitation  inception  is  very  sensitive  to  Reynolds  scaling.  McCormick’s  experimental 
data  [1954]  show  that  cavitation  inception  is  well  represented  by 


aif/aim=  (Ref/Rem)n.  (7) 

Shen,  Jessup  and  Gowing  [2001]  have  derived  a  theory  for  high  Reynolds  number  application. 
Their  theoretical  result  gives 

5.16  Log  (Log  Ref/LogRem)/ 

/  Log  (Ref/Rem)  ^ 

The  numerical  values  of  n  as  function  of  Reynolds  numbers  are  given  in  Figures  3  and  4.  The 
values  of  n  decrease  as  the  values  of  Reynolds  number  increase.  For  examples,  n  =  0.40,  0.33, 
and  0.31  at  Re  =  3x1 0s,  5xl06,  and  107,  respectively. 

Based  on  similarity  flow  approaches,  Shen,  Jessup  and  Gowing  [2001]  have  theoretically 
derived  the  ratio  of  vortex  core  sizes  in  model  and  full-scale  as 

acf/acm  =  X(Ref/Rem)n/2,  (9) 

where  X  is  the  linear  scale  ratio  of  full-scale  and  model,  Rem  and  Rerare  model  and  full-scale 
Reynolds  number,  respectively.  Equation  (9)  shows  the  viscous  effects  on  the  vortex  core  size. 
As  to  be  shown  later  in  the  numerical  section,  the  ratio  of  maximum  bubble  radii  computed  by 
Hsiao  et  al  [2000]  are  related  to  the  vortex  core  sizes  as  following 

Rjriax  m  /  acm  —  Rmax  f  /  .  (10) 
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In  dimensional  analysis  for  scaling,  a  variable  must  be  expressed  in  a  non-dimensional 
form.  As  a  cavitation  bubble  is  convected  downstream  inside  the  trailing  vortex  core,  the 
pressure  field  in  the  vortex  determines  the  growth  and  collapse  of  the  bubble.  For  the  pressure 
distributions  inside  a  vortex  core  to  be  represented  by  a  Rankine  model,  Equation  (10)  states  that 
the  characteristic  length  to  normalize  the  maximum  bubble  radius  is  the  vortex  core  size.  From 
Equations  (9)  and  ( 1 0),  we  have 

Rmaxm/Rmaxf  =(Ref/Rem)n/2/k.  (11) 

The  ratio  of  dynamic  variables  Rmax  m  /  Rmax  f  has  now  been  expressed  in  geometric  and 
kinematic  variables. 


The  Ratio  of  ^  e "  ^  f  c  "P(  ))m  jn  a  jrajjjng  Vortex:  The  pressure  in  a 
trailing  vortex  for  a  bubble  to  collapse  is  defined  in  Equation  (1).  At  R=  Rmax,  the  terms  pg  and 
2S/Rmax  are  small.  Recall  that  Pe  is  the  encounter  pressure  inside  the  vortex  core  local  to  the 
bubble  at  R  =  Rmax.  For  a  vaporous  cavitation  bubble,  Equation  (1)  gives: 

p(R)-  Pe  *  pv  -  Pe  =  pv  -  Poo  +  Poo  -  Pe 

=  -  (0.5pVoc2)  (Cp  +  cr ) ,  (12) 

where  Cp s (Pe - P* )/ (0.5 pV^2)  and  as(P,- pv )/ (O.SpV^2)  .  (13) 

So 

V(P<R)-r.)r  iJWhfX  -  f(c,  +4  /V-(c,  +o).(vrt/v„) 

where  pf/pm  has  been  set  to  equal  1 . 

Assume  that  cavitation  inception  occurs  when  the  pressure  in  the  core  drops  below  the 
vapor  pressure.  We  have: 


C7jf  /  Ojiyi  Cpminf  /  Cpminin  • 

From  a  Rankine  vortex  model,  Equations  (12),  (13)  and  (15)  give 

V(P(R)-Pe)f  /V(p(R)-Pe)m  =(Ref/Rem)^(Vo,f/V00m) 

The  ratio  of  dynamic  variables 

V(P(R)  Pc)f/V(P(R)-Pe)m 

at  R  =  Rmax  has  now  been  expressed  in  geometric  and  kinematic  variables. 


(15) 


(16) 
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SCALING  OF  FREQUENCY  AT  PEAK  AMPLITUDE 


For  pf/pm  =  1 ,  Equation  (6)  gives 


f  =f  fR.,..YV<p(R)-'>.), 


m  J 


=  f  - 
pm  x 


Ref 

Re! 


\V 


V. 


xf 


vY™. 


(17) 


As  shown  in  Equation  (8),  at  very  high  Reynolds  numbers,  the  value  of  n  approaches 
zero.  Only  in  this  condition,  Equation  (17)  gives  the  same  result  as  Strasberg’s  well-known  X 
frequency  scaling  [1957,  1977]  if  viscous  flow  effects  are  neglected.  Namely,  Equation  (17) 
provides  tip  vortex  Reynolds  number  correction  to  Strasberg’s  X  frequency  scaling. 

Effect  of  Free  Stream  Nuclei  Size  on  Frequency  Scaling :  As  shown  in  Equation  (1), 
when  a  bubble  is  in  a  static  equilibrium,  the  balance  of  pressure  across  the  bubble  interface  can 
be  written: 


Pe  =  Pv+Pg-2S/R.  (18) 

The  bubble  becomes  unstable  and  experiences  rapid  growth  when  the  derivative  of  pressure 
gradient,  dP/dR  >  0.  Because  of  the  intrinsic  property  of  bubble  dynamics,  at  a  given  cavitation 
number,  a  bubble  grows  to  about  the  same  maximum  size  regardless  of  the  initial  bubble  sizes 
(Chahine  and  Shen  [1986],  Strasberg  [1957]).  The  frequency  scaling  given  in  Equation  (17) 
should  not  be  affected  by  free  stream  nuclei  size  distributions  as  long  as  the  cavitation  noise  is 
generated  by  non-interacting  bubbles. 


SCALING  OF  ACOUSTIC  AMPLITUDE 

Acoustic  pressure :  The  relation  between  the  radial  motion  of  a  growing  and  collapsing 
cavity  and  the  pressure  in  the  surrounding  liquid  have  been  treated  extensively  in  the  literature 
relating  to  cavitation  and  underwater  explosion  (see  Strasberg  [1956],  Fitzpatrick  [1958]).  For 
incompressible  flow  we  have  (Fitzpatrick  [1958]): 

p(r,t)  — P e(t)  =  [  pR2  d2R  /dt2  +  2pR(dR/dt)2]/r  -  [  0.5pR4(dR/dt)2/r4]  .  (19) 

where  p(r,t)  is  the  pressure  in  the  liquid  at  radial  coordinate  r  and  time  t;  R  is  the  bubble  radius; 
and  Pe(t)  is  the  pressure  inside  the  vortex  core  local  to  the  bubble  at  time  t.  It  is  noted  that  on  the 
bubble  surface  where  r  =  R,  Equation  (19)  reduces  to  Equation  (1)  when  p(R,t)  denotes  the 
pressure  in  the  liquid  just  outside  the  bubble  surface. 
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At  a  sufficiently  great  distance  from  the  cavity,  r»R,  the  second  term  on  the  right  hand 
side  of  Equation  (19)  is  negligible.  The  acoustic  pressure  Pa  =  p(r,t)  -  Pe  (t)  can  be  expressed  by 

Pa  =  pR  [R  d2R  /dt 2  +  2(dR/dt)2]/r  =  p/47tr[d2o  /dt2] ,  (20) 

where  o  =47tR3/3,  denotes  the  volume  of  the  bubble.  Equation  (20)  is  the  well-known  acoustic 
pressure  formula,  which  states  that  the  acoustic  pressure  is  related  to  the  second  time  derivative 
of  volume  change  and  the  inverse  of  the  distance  from  the  noise  source. 

Peak  Acoustic  Amplitude :  A  sketch  of  an  acoustic  spectrum  is  shown  in  Figure  2.  It  is 
well  known  (Knapp  et  al  [1970])  that  the  peak  amplitude  of  noise  spectrum  is  related  to  the 
maximum  bubble  radius,  R™*  .  Let  Papeak  denotes  the  peak  acoustic  amplitude  at  R  =  Rmax.  We 
have 

dR/dt  =  0  at  R  =  R  max  (21) 

Substitute  Equation  (21)  into  Equation  (20),  we  obtain 

Papeak  =  pRmax2d2R/dt2/r.  (22) 

At  R  =  Rmax ,  dR/dt  =  0  and  Equation  (19)  gives 

d2R  /dt2  =  (1/Rmax)  (pv  +  pg-  Pe  -  2S/R)/  p  .  (23) 

From  Equations  (22)  and  (23),  the  peak  acoustic  amplitude  can  be  expressed  by 
Papeak  (  Rmax  /  f  )  (Pv  4"  Pg —  Pe(Rmax)  2S/Rmax) 


=  (R  max  /r)  (p(R  max  )-Pe(R  max))  • 


(24) 


If  the  shapes  of  the  acoustic  spectra  due  to  bubble  cavitation  are  similar  between  the 
model  and  full-scale,  scaling  of  the  acoustic  amplitude  can  be  expressed  by 


Paf  /  Pam  Papeak  f  /  Papeak  m 

=  (R  max  f/Rmax  m  )(^m  I  ff)  [p(R  max)  Pe(Rmax)]f/[p(Rmax)-Pe(Rmax)]m  •  (25) 

Another  approach  to  derive  Equation  (25)  is  by  dimensional  analysis  as  discussed  by 
Strasberg  [1977].  We  know  that  the  intensity  of  acoustic  pressure  is  inversely  proportional  to  the 
distance  from  the  noise  source.  We  need  a  characteristic  length  scale  to  normalize  the  distance: 
that  length  scale  is  the  maximum  radius.  We  need  a  characteristic  pressure  to  normalize  the 
acoustic  pressure:  that  pressure  scale  is  the  driving  pressure,  p(R)  -  Pe,  responsible  for  bubble 
growth  and  collapse. 
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Recall  that 


Rmaxm/Rmaxf  =  (Ref/Rem)"/2/X.  (11) 

and 

[p(Rmax)  -  Pe(Rmax)]f /[p(Rmax)  -  Pe(Rmax)]m  =  (Ref  /  Rem  )"  (VKf/VMm)2  .  (16) 

From  Equations  (1 1),  (16)  and  (25),  we  obtain  the  amplitude  scaling  as  follows  : 

Par  /  Pam  =  [X  (Ref  /  Rem  )'n/2]  (rm  /  rf)  [Ref  /  Rem  )n  (V«f . /VMm)2] 

=  X  (Ref  /  Rem  f2  (rm  /  rf)  (Vxf/Vxm)2  .  (26) 

The  analysis  so  far  has  assumed  that  similar  frequency  bands  will  be  used  for  both  model  and 
prototype.  Following  Strasberg’s  approach,  the  relationship  between  noise  levels  Lf  and  Lm  for 
prototype  and  model  is 

Lf  =  Lm  +  20  Log  X  +  10  Log  (Ref/Rem  )n  +  20  Log  (rm  /rf)  +  20  Log  (V«f /V^)2  •  (27) 

The  analysis  so  far  is  related  to  a  single  cavitation  bubble.  In  case  there  are  multiple 
events  of  bubble  cavitation,  Equation  (27)  can  be  modified  to  include  multiple  events  as  long  as 
the  events  are  incoherent  and  the  pressures  from  all  signals  add  up.  Let  N  denote  the  number  of 
cavitation  events  per  unit  time.  We  have 

Lf  =  Lm  +  20  Log  X  +  1 0  Log  (Ref  /Rem  )n  +  20  Log  (rm  /rf) 

+  20  Log  (V«,f /VTOtn)2  +  1 0  Log  (Nf/Nm)  .  (28) 

Comparison  with  Strasberg  Noise  Level  Scaling :  Using  dimensional  analysis,  Strasberg 
has  derived  a  noise  level  scaling  as  follow.  Let  D  denote  a  length  scale  such  as  a  propeller 
diameter,  and  P«  the  ambient  pressure.  He  has  obtained  the  non-dimensional  n  parameter, 

re  =  (Pa/Poo)(r/D) 


or 

Paf  /  Pam  —  (Df  /  Dm  )(rm  /  rf )  (Poof  /  Poom  )  •  (29) 

A  comparison  of  Equations  (25)  and  (29)  shows  that  amplitude  scaling  given  in 
Equation  (29)  by  Strasberg  and  Equation  (26)  by  the  present  derivation  will  be  the  identical  if 


(Df  /  Dm  )  Rmax 


f/R 


max  m  ♦ 


(30a) 


and 


(Poof  /  Poon,  )  =  (p(R)  -  Pe)f  /  (P(R)  -  Pe)m  - 


(30b) 
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Equations  (30a)  and  (30b)  show  that  in  the  dimensional  analysis,  if  the  maximum  bubble  radius 
is  used  as  the  characteristic  length;  and  the  driving  pressure  inside  the  vortex  core  is  used  to 
normalize  the  acoustic  pressure,  then  the  present  scaling  and  Strasberg’s  scaling  on  acoustic 
amplitude  are  identical. 


CAVITATION  NUMBER  SCALING 

The  present  theoretical  derivation  is  based  on  a  similarity  flow  approach.  This  means 
that  pressure  fields  and  cavitation  patterns  in  vortex  cores  are  similar  between  full-scale  and 
model.  This  similarity  requirement  is  satisfied  by  setting 


Of  /  =  (Patm  P§h  —  Pv)f /(Patm  P§h  pv)m  (Ref  /  Rem  )  ,  (31) 

where  Patm  denotes  the  atmospheric  pressure  on  the  water  surface  and  h  denotes  the  depth  of 
submergence.  Since,  pv  «  Patm  and  pf  ~  pm  =  p,  we  have 

hf  =  [(Patm  +  pghm  )  (Ref  /  Rem  )n  -  Patm)]  /  pg .  (32) 

Equation  (32)  states  that  frequency  and  amplitude  scaling  of  acoustic  spectrum  measured  at  a 
water  depth,  hm,  are  to  be  applied  to  full-scale  by  Equations  (17)  and  (28)  at  the  full-scale  depth 
hf  given  by  Equation  (32). 


NUMERICAL  VALIDATION  OF  THE  ACOUSTIC  SPECTRA  SCALING 

Two  geometrically  similar  propellers  of  different  sizes  are  considered  [see  Figure  1]. 
The  two  propellers  have  Rp  =  61  cm  and  244  cm,  respectively.  The  small  propeller  is  termed 
model  propeller  and  the  large  propeller  is  termed  full-scale  propeller  in  this  numerical  study. 
The  geometric  scale  ratio  X  of  these  two  propellers  is  X,  =  4.  In  the  present  numerical  study,  the 
upstream  velocity  is  assumed  to  be  the  same  for  both  full-scale  and  model,  Voof  =  V«,m.  The  axial 
velocity  in  the  vortex  core  is  assumed  to  be  13  m/s  for  both  full-scale  and  model.  The 
assumption  of  equal  axial  velocity  between  full-scale  and  model  is  reasonable. 

Vortex  core  size  ac,  circulation  T  and  tangential  velocity  distributions  Vt  along  the  x-axis 
on  a  trailing  vortex  has  been  measured  with  LDV  in  36-inch  water  tunnel  on  an  open  propeller 
(Judge  et  al  [2001]).  Using  the  scaling  formula  discussed  in  a  paper  by  Shen  et  al  [2001],  the  36- 
inch  water  tunnel  data  are  scaled  to  the  model  propeller  used  in  this  numerical  study  and  given  in 
Figure  5.  The  pressure  distributions  are  computed  by  using  the  Rankine  vortex  model.  The 
hydrodynamic  data  on  the  full-scale  propeller  are  given  in  Figure  6  by  applying  McCormick’s 
hypothesis 


/  cm  -  Cpf  /  -  Cpm  (Ref  /  Rem  )  .  (33) 

In  the  present  numerical  study,  we  use  n  =  0.4  and  the  same  kinematic  viscosity  of  water  for 
model  and  full-scale. 
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Acoustic  spectrums  due  to  bubble  growth  and  collapse  in  a  trailing  vortex  core  are 
numerically  computed  to  compare  with  the  scaling  formula  derived  from  a  similarity  flow 
approaches  as  shown  in  the  above  equations. 


Bubble  Dynamic  Model 

It  is  assumed  that  the  bubble  entering  the  vortex  core  remains  spherical  throughout  the 
cavitation  process.  Bubble  dynamics  in  the  cavitation  process  can  be  computed  by  the  SAP 
modified  Rayleigh-Plesset  Equation  (Hsiao  et  al  [2000]) 

R  (d2R/dt2)  +  3/2  (dR/dt)2  =  1  /p  (pv  +  pg  -  Pe  -  2S/R  -  4p(dR/dt)/R)  +  p(0  -  Ub)2/4),  (34) 

where  0  -  0b  denotes  the  slip  velocity  between  the  fluid  and  the  bubble.  The  viscosity  effect  in 
Rayleigh-Plesset  Equation  is  also  included  in  Equation  (34),  which  is  used  in  the  present 
numerical  calculations.  By  taking  the  Fourier  transform  of  the  radiated  acoustic  pressure  Pa 
given  in  Equation  (20),  the  acoustic  spectra  can  be  obtained. 

Numerical  results:  A  small  portion  of  numerical  computations  is  provided  in  this  report 
for  discussions.  The  history  of  bubble  growth  and  collapse  for  the  model  and  full-scale 
conditions  are  given  in  Figures  7  and  8  for  the  initial  bubble  size  of  Ro  =  5  pm.  Consider  first 
the  model  case  of  Figure  7.  At  om  =  7.50  (upper  figure),  the  5  pm  bubble  only  experiences 
marginal  growth  as  the  bubble  is  convected  downstream  with  very  small  change  in  acoustic 
pressure.  As  the  ambient  pressure  is  slightly  reduced  to  am  =  7.49  (lower  figure),  the  bubble 
grows  to  a  maximum  size  of  Rmax  =195  pm  at  time  t  around  0.0081  sec.  The  bubble  then  enters  a 
high-  pressure  zone  and  the  bubble  begins  to  collapse  and  rebound.  The  acoustic  pressure 
initially  rises  above  the  background  level  of  zero  then  drops  to  negative  as  the  bubble  grows  to  a 
maximum  size.  The  acoustic  pressure  rises  to  positive  as  the  bubble  begins  to  collapse.  A  sharp 
rise  in  acoustic  pressure  is  computed  when  the  bubble  size  reaches  the  minimum.  Since 
compressibility  due  to  the  presence  of  un-dissolved  gas  inside  the  bubble  is  included  in  this 
bubble  dynamic  model,  the  minimum  bubble  size  at  the  collapse  is  finite.  This  follows  growth 
and  rebound  as  shown  in  Figure  7.  The  qualitative  description  of  rise  and  fall  of  acoustic 
pressure  in  response  to  bubble  growth  and  collapse  is  in  agreement  with  the  results  given  in 
Fitzpatrick  and  Strasberg  paper  [1958].  Note  that  at  time  t  =  0,  the  bubble  is  located  at  the 
propeller  blade  tip.  The  history  of  pressure  Pe,  in  the  vortex  core  to  be  encountered  by  the 
bubble  and  the  change  in  gas  pressure  pg  inside  the  bubble  are  also  shown  in  Figures  7  and  8. 

Next  consider  the  full-scale  case  of  Figure  8  for  the  initial  bubble  size  of  5  pm.  Due  to 
Reynolds  scale  effect  on  vortex  core  pressure  distributions,  cavitation  inception  occur  at  Of  = 
13.09.  The  rise  and  fall  of  acoustic  pressure  in  response  to  bubble  growth  and  collapse  along  the 
trailing  vortex  show  the  same  patterns  in  the  full-scale  and  in  the  model. 


Scaling  of  Maximum  Bubble  Radius 

Numerical  Results:  The  maximum  bubble  radius  Rmax  vs  normalized  cavitation  number  a 
is  shown  in  Figures  9  and  10  for  the  initial  bubble  radii  of  5  and  10  pm,  respectively.  The 
cavitation  numbers  are  normalized  by  (Rem  /  Re)0  4,  where  Re  =  Rem  for  the  model  scale  curve 
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and  Re  =  Ref  for  the  full-scale  curve.  Cavitation  inception  occurs  around  the  normalized 
cavitation  numbers  of  7.50  and  7.52.  Once  the  cavitation  is  beyond  inception  (o  <  7.5),  the 
maximum  bubble  size  increase  sharply  and  almost  linearly  with  the  reduction  of  cavitation 
number.  Instead  of  Rmax,  we  can  plot  Rmax  f  /  Rmax  m  vs  normalized  a  as  shown  in  Figure  11. 
Except  in  the  neighborhood  of  cavitation  inception  region,  the  bubble  dynamic  calculations  give 
the  ratio  of 


Rmax  f  /  Rirrax  m  ~  3.0  .  (35a) 

Present  Similarity  flow  method:  The  ratio  of  maximum  bubble  radius  derived  from 
similarity  flow  approach  is  given  in  Equations  (10)  and  (11): 

Rmax  f /  Rmax  rn  =  acf  /  acm  =  l  (Ref  /  Rem  f2  =  4°  8  =  3.03  .  (35b) 

This  value  of  3.03  derived  from  similarity  flow  approach  is  in  very  good  agreement  with  the 
value  3.0  obtained  from  bubble  dynamics  numerical  calculations.  The  good  agreement  between 
similarity  flow  method  [Equation  (35b)]  and  numerical  calculations  [Equation  (35a)]  supports 
the  assumption  used  in  dimensional  analysis  that  the  characteristic  length  to  normalize  the 
maximum  bubble  size  is  the  vortex  core  radius  [Equation  (10)],  not  a  propeller  diameter  or  a 
blade  chord  length. 


Scaling  of  Frequency  at  Peak  Amplitude  of  Acoustic  Spectrum 

The  acoustic  pressure  is  computed  from  Equation  (20).  The  acoustic  amplitude  spectrum 
is  computed  from  the  following  equation. 


s  = 


Pa  e‘i,otdt 


12  n 


(36) 


Acoustic  spectra  have  been  calculated  for  several  initial  bubble  radii  and  cavitation 
numbers.  As  an  example,  the  spectra  expressed  in  a  narrow  band  analysis  are  shown  in  Figure 
12  with  the  initial  bubble  radius  of  10  pm  at  a  =  7.2  for  the  model  and  a  =  7.2  x  (4)0  4  =  12.53  for 
the  full-scale,  respectively.  In  the  plots,  the  notation  D  =  0.5  ac  denotes  that  the  1 0  pm  bubble  is 
released  at  half  radial  distance  from  vortex  core  center  to  the  edge  of  the  vortex  core.  The 
bubble  is  released  at  the  upstream  location  of  x  =  0,  corresponding  to  the  blade  tip  location. 
Discussion  of  the  dynamic  process  of  a  bubble  spiraled  into  the  vortex  center  can  be  found  in  the 
report  by  Hsiao  et  al  [2000].  Another  example  is  given  in  Figure  13  for  the  initial  bubble  radius 
of  2.5  pm  at  a  =  7. 1 5  for  the  model  and  o  =  7. 1 5  x  (4)0'4  =  12.45  for  the  full-scale,  respectively. 

Numerical  results :  Figures  12  and  13  give  the  frequency  at  peak  amplitude  to  be  around 
fPf  ~  900  Hz  for  the  full-scale  spectrum  and  fpm  ~  2 1 1 0  Hz  for  the  model  spectrum  for  both  10 
and  2.5  pm  bubbles.  Both  full-scale  and  model  frequencies  are  relatively  low  when  compared 
with  some  of  the  water  tunnel  measurements.  This  is  because,  we  only  consider  cases  where  the 
maximum  pressure  peak  occurs  at  Rmax.  It  is  remarked  that  the  peak  frequency  depends  on  the 
vortex  core  size  and  the  encounter  pressure.  The  peak  frequency  can  be  changed  if  these  two 
parameters  are  changed.  Numerical  calculations  give 
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fpf  -  fpm  (900  /2110)  =  fpm  /2.34  . 


(37a) 


This  numerical  result  is  quite  different  from  Strasberg’s  classic  X  scaling.  Note  that  the  linear 
scale  ratio  used  in  this  study  is  X  =  4. 


Present  Similarity  flow  method:  Equation  ( 1 7)  gives  the  scaling  formula  by 


f  =f  I  — 

*pf  Apml  * 


Ref 


YY 


Re 


m  J 


V. 


V. 


Note  that  V«,f  =  ,  Vf  =  vm,  X  =  4,  and  n  =  0.4  are  used  in  this  numerical  study. 


(37b) 


The  good  agreement  between  the  similarity  flow  method  [Equation  (37b)]  and  the 
numerical  calculations  [Equation  (37a)]  supports  the  assumption  used  in  dimensional  analysis 
that  the  characteristic  length  and  pressure  term  to  normalize  the  peak  frequency  are  maximum 
bubble  radius  and  encounter  pressure  in  the  vortex  core.  All  other  calculations  not  shown  in  this 
report  give  the  same  good  agreement  between  frequency  scaling  of  Equation  (17)  and  bubble 
dynamic  calculations. 


Numerical  calculations  show  that  for  a  given  cavitation  number,  all  bubbles  of  different 
initial  radii  grow  to  about  the  same  maximum  bubble  radius  Rmax.  This  is  an  intrinsic  property  of 
the  bubble  dynamic  equation.  This  means  that  frequency-scaling  formula  given  in  Equation  (17) 
holds  regardless  of  the  initial  bubble  size  in  the  free  stream. 


Scaling  of  Peak  Amplitude  of  Acoustic  Spectrum 

Numerical  Results :  Figure  12  gives  the  model  and  full-scale  amplitude  spectra  for  the 
initial  bubble  radius  of  Ro  =10  pm.  The  computed  peak  amplitudes  for  the  full-scale  and  model 
spectra  are  pf  a  2.5  x  10‘2  Pa  and  pm  a  4.5  x  10“3  Pa  for  the  model,  respectively.  The  distances 
between  the  hydrophone  location  and  the  cavitation  bubble  used  in  these  calculations  are  rr  = 
0.52  m  for  the  full-scale  and  rm  =  0.32  m  for  the  model. 

Figure  13  gives  the  model  and  full-scale  amplitude  spectra  for  the  initial  bubble  radius  of 
Ro  =  2.5  pm.  The  peak  amplitudes  for  the  full-scale  and  model  spectra  are  pf  a  7.0  x  10‘2  Pa  and 
pm  «  1.4  x  10'2  Pa  for  the  model,  respectively.  The  distances  between  the  hydrophone  location 
and  the  cavitation  bubble  are  rf  =  0.52  m  for  the  full-scale  and  rm  =  0.32  m  for  the  model. 

Let  Lf  and  Lm  denotes  the  noise  level  generated  by  the  cavitation  bubble  for  the  full-scale 
and  model,  respectively.  The  differences  in  noise  level  between  full-scale  and  model  for  Ro  =  10 
and  2.5  pm  are  given  by 

AL  =  Lf  -  Lm  =  20  Log  (pf  /  pm) 

=  20  Log  (2.5  x  10'2  /  4.5  x  10'3  )=  14.9  dB  forR0=10pm  (38a) 

=  20  Log  (7.0  x  10'2/ 1.4  x  10’2)=  14.0  dB  forR0  =  2.5pm.  (38b) 

12 


Present  Similarity  flow  method:  The  amplitude  scaling  is  given  in  Equation  (28)  as 
follows 

AL  =  Lf  -  Lm  =  20  Log  X +  1 0  Log  (Ref  /Rem  )n  +  20  Log  (rm  /rf) 

+  20  Log  (Vccf  AU,)2  +  10  Log  (Nf/Nm) 

=  20  Log  4  +  10  Log  (4) 04  +  20  Log  (0.32  /  0.52)  =  10.3  dB.  (38c) 

Note  Nf  =  Nm  =  1  for  a  single  bubble  study. 

A  comparison  of  numerical  calculations  and  amplitude  scaling  formula  given  in 
Equations  (38a),  (38b)  and  (38c)  indicates  a  difference  of  3.7  dB  for  R0  =  2.5  pm  and  4.6  dB  for 
R0  =  10  pm.  Analysis  of  the  numerical  data  indicates  that  the  cavitation  bubble  experiences 
growth  and  collapse  at  slightly  different  normalized  encounter  pressure  between  the  model  and 
full-scale.  This  slight  difference  in  encounter  pressure  has  very  minor  effect  on  frequency  but 
some  effect  on  acoustic  amplitude.  Further  study  is  needed  to  quantify  the  initial  bubble  size 
effect  on  acoustic  amplitude.  At  present.  Equation  (28)  is  considered  to  be  valid  in  a  global 
sense  with  uncertainty  of  a  few  dB. 


CONCLUSIONS 

Scaling  of  noise  generation  due  to  bubble  cavitation  in  a  trailing  vortex  has  been 
theoretically  investigated.  Three  scaling  formulas  have  been  derived  based  on  geometric, 
kinematic  and  dynamic  similarities  between  full-scale  and  model.  The  first  two  scaling  formulas 
provide  frequency  and  noise  level  (amplitude)  scaling  to  obtain  full-scale  spectrum  from  model 
measurements.  The  third  cavitation  number  scaling  relates  the  depths  or  speeds  between  full- 
scale  and  model  for  application  of  the  acoustic  scaling. 

Due  to  Reynolds  scaling  effects  on  vortex  core  dynamics  and  pressure  distributions  in  a 
trailing  vortex,  the  present  frequency  scaling  is  noticeably  different  from  Strasberg’s  classic  X 
scaling.  However,  if  viscous  effects  are  not  considered,  the  present  theory  gives  the  same 
frequency  scaling  as  Strasberg’s  X  scaling.  The  present  frequency  scaling  represents  a  Reynolds 
scaling  correction  of  Strasberg’s  X  scaling. 

The  frequency  scaling  is  found  to  be  insensitive  to  nuclei  sizes. 

The  present  noise  level  (amplitude)  scaling  also  represents  a  Reynolds  scaling  correction 
of  Strasberg’s  noise  level  scaling.  If  viscous  effects  are  not  considered,  the  present  noise  level 
scaling  agrees  with  Strasberg’s  noise  level  scaling. 

Acoustic  narrow  band  spectrums  generated  by  a  spherical  cavitation  bubble  are 
numerically  studied  by  using  a  bubble  dynamic  computer  code  developed  by  Dynaflow  Inc. 
Frequency  scaling  formula  derived  from  similarity  flow  approach  is  shown  to  be  in  excellent 
agreement  with  bubble  dynamic  calculations. 
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Noise  level  (amplitude)  scaling  formula  derived  from  similarity  flow  approach  is  found 
to  differ  3.7  dB  for  Ro  =  2.5  pm  and  4.6  dB  for  Ro  =  10  pm  from  bubble  dynamic  calculations. 

Further  numerical  studies  are  needed  to  quantify  nuclei  size  effect  on  noise  level  scaling. 
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Cavitation  Spectrum  (Ross  1976) 
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Figure  3.  The  Values  of  n  Computed  vs.  Measured 
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Figure  4.  The  Values  of  n  for  1/4-scale  and  1/16-scale  Models 


Figure  5.  Pressure  Distributions  in'the'Vortex  Core  used 
in  the  Present  Numerical  Calculations  :  Model 


Figure  6.  Pressure  Distributions  in  the  Vortex  Core  used 

in  the  Present  Numerical  Calculations  :  Full-Scale 
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Figure  7.  Computed  Bubble  Radius  History,  Gas  Pressure,  Encountered  Pressure  and 
Radiated  Acoustic  Pressure  at  Ro  =  5  pm  and  a  =  7.50  and  7.49  for  Model  Scale 
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Figure  8.  Computed  Bubble  Radius  History,  Gas  Pressure,  Encountered  Pressure  and 
Radiated  Acoustic  Pressure  at  Ro  =  5  pm  and  <7=13.10  and  13.09  for  Full-Scale 
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Figure  9.  Maximum  Bubble  Size  versus  Normalized  Cavitation  Number 
for  Full-Scale  and  Model  with  Ro  =  5  (am 


Figure  10.  Maximum  Bubble  Size  versus  Normalized  Cavitation  Number 
for  Full-Scale  and  Model  with  Ro  =  10  jam 
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Figure  1 1 .  R  max  f  /  Rmax  m  VS.  (7 


Figure  12.  Numerically  Computed  Acoustic  Spectra  for 
Initial  Bubble  Size  of  Ro  =  1 0  pm 


23 


Frequency (hz) 

Figure  13.  Numerically  Computed  Acoustic  Spectra  for 
Initial  Bubble  Size  of  Ro  =  2.5  pm 
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